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Abstract 

We study asymptotics of fiber integrals depending on a large parameter. When the 
critical fiber is singular, full-asymptotic expansions are established in two different 
cases : local extremum and isolated real principal type singularities. The main co¬ 
efficients are computed and invariantly expressed. In the most singular cases it is 
shown that the leading term of the expansion is related to invariant measures on the 
spherical blow-up of the singularity. The results can be applied to certain degenerate 
oscillatory integrals which occur in spectral analysis and quantum mechanics. 
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1 Introduction and statement of the main result. 

In [3] J.Bruning& R.Seeley have studied asymptotic expansions of integrals : 

OO 

H{z) = j a{xz, x)dx, z ^ oo, (1) 

0 

where a{x, is a singular symbol. The result of [3] is quite remarkable, in par¬ 
ticular because it can be directly applied to spectral analysis. Many asymptotic 
questions can be reduced to the study of the previous problem but it is also 
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interesting to consider a generalization : 


I{z) = J g{zf{x),x)dx, z 

X 


oo. 


( 2 ) 


where g : Wx X ^ W, f : X ^ M are smooth and X is a smooth differentiable 
manifold equipped with the C°° density dx. As a motivation, we remark that 
the trace formula for certain linear self-adjoint operators can be stated in the 
form of Eq. (2). In the semi-classical regime when / is a quadratic form this 
problem was involved in [2], See also [5] for the case of / with an isolated and 
degenerate singularity associated to an homogeneous dehnite jet. 

General assumptions. Throughout this work, we will assume that |/| is strictly 
positive outside of a compact set and that the Fourier transform g w.r.t. t exists 
with d^g{t,x) G x X), \/k. For this reason g will be called a symbol. □ 

Remark 1 This assumption on g is strong but can be weakened. More general 
conditions can be found in [7]. Mainly, this condition will be used to reach 
integrals with compact supports. 

As cx), the asymptotic behavior of I{z) is related to the critical fiber : 

6 = r\m = {xex / fix) = 0 }. ( 3 ) 


This can easily be viewed with the Fourier inversion formula : 

I{z) = J J e^^^^'^^^git,x)dtdx, z ^ oo, 

X R 

where g{t,x) is the normalized Fourier transform of g w.r.t. t : 

e"*^*c/(r,a;)dr, 

R 



( 4 ) 


( 5 ) 


In Eq. (4), the stationary points w.r.t. t are precisely given by ©, i.e. liz) is 
asymptotically supported by ©. Since / is smooth © is closed and according 
to the general assumptions above we obtain : 

(i7o) The fiber & is compact. (6) 

(Rq) is not absolutely necessary, e.g. if g decreases fast enough near the bound¬ 
ary of X n G. But, to simplify, we will only consider the compact case. To 
obtain an easier formulation of the problem we recall an elementary result. 

Lemma 2 If d^g G x X), V/c G N, modulo terms Oiz~°^), asymptotics 
of Eq. (2) are not changed by assuming that g is compactly supported near &. 
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Proof. With © compact we choose a cut-off function 0 G C“(X) such that 
\1/ = 1 near © and 0 < \1/ < 1. We shall estimate the error integral : 

E{z) = J J —'^{x))dtdx, z ^ oo. (7) 

X R 

With L = —{i/zf{x))dt, we have V/c G N. By integration 

by parts and since |/(a;)| >C on supp(l — d'), we obtain : 

\Eiz)\ < {Cz)-^\\d^,g{t,^)\W[^>.x) = 0{z-^), \/k G N. (8) 

This gives the desired result, with our hypothesis on g. ■ 

Lemma 1 allows to consider only integrals with compact support w.r.t. x 
which simplihes all questions of convergence. Notice that we can weaken the 
condition on g to d^g G L^(M x X), Wk < k^ with an error 0(z~’^°). We 
are mainly interested in the sitnation where © has an isolated singularity. 
If xo G © is such a critical point, let 0 G C^{X) be a cnt-off microlocally 
supported near xq. We split-up our integral as I{z) = Ir{z) + Is{z), where : 

^r{z) = J g{zf{x),x){l - e){x)dx, (9) 

Is{z) = J g{zf{x)t,x)Q{x)dx. (10) 

X 

This procedure can be extended with hnitely many critical points on ©. The 
regular part 1^ can be treated by the generalized stationary phase method, 
with non-degenerate normal Hessian, which we recall in section 2. Since Ig is a 
local object and the main contribntions below concern invariant objects, there 
is no loss of generality to assume that supp( 0 ) is an open of M”, n = dim(X). 
For xq G © a singularity of hnite order we can write the germ of / as : 

/(i)=A(i)+o(ii(i-i„)ir‘), (11) 

where 7 ^ 0 is homogeneous of degree k > 2 w.r.t. {x — xq). The first 
elementary result concerns extremum attached to such homogeneous germs : 

Theorem S If f has a local extremum xq on © whose jet is given by Eg. (11) 
(a fortiori k is even), we obtain a full-asymptotic expansion : 

(12) 

JGN 

//dim(X) = n, the leading term is given by : 

Is{z) = z~'^ (tl^,g{t,xo)'^^ J \fk{9)\~'^d9 + 0{z~^), (13) 

§n-l 
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with te = max(t, 0) if xq is a minimum and max(—1,0) for a maximum. 

The reader can observe that Theorem 3 includes the case of a non-integrable 
singularity on © for k > n. Accordingly, always for k > n, the contribution of 
the critical point is bigger than the regular contribution (see section 2). 

Remark 4 The result stated in Eq. (13) includes the case k = 2 and the 
integral on the sphere can be expressed in terms of special functions. Under 
ad-hoc conditions, Theorem 3 can be extended to a singularity which is a sum 
of positively homogeneous singularities, e.g. f(xi,X 2 ) = ||a^i||"^ + ||a^ 2 ||®- 

The case of non-extremum degenerate critical points is more difficult. Since 
this problem can be very complicated in general position we impose : 

{Hi) & has a unique critical point xq. Moreover, fk defined in Eq. (11) is 
non-degenerate in the sense that : 

Vfk ^ 0 on n {h = 0}. (14) 

By homogeneity, {Hi) is equivalent to ” fk has an isolated singularity”. Observe 
that this condition is very close to Hormander’s real principal type condition 
for distributions. We dehne the integrated density of fk on the sphere (or 
co-area) as : 


LVoKw) = J \dL\{e), dL{e) Adfk{9) = de. (15) 

{fk(S)=w} 


Where dL is the n—2 dimensional Liouville measure induced by fk on i.e. 
the Riemannian density induced by fk on the standard density of Note 
that {Hi) insures that LVol(tc) is well dehned and smooth near the origin. 

Theorem 5 Under the previous assumptions and if xq satisfies {Hi), the sin¬ 
gular part of our integral admits a full asymptotic expansion : 

OO . OO 

^s{z) ~ log(^). (16) 

i=o j=o 

a) If k > n (non-integrable singularity), the leading term is : 

h{z)=Coz-"^+0{z-^\og{z)), (17) 

where the distributional coefficient Cq is given by : 

^{(tl~^,g{t,xo)'^ J \fk{0)\~'^de + (tI~^,g{t,xo)'^ J \fk{9)\-'^ dO). 

{/fc>0} {/fc<0} 
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b) If k divides n the leading term is logarithmic : 


Is{z) = Dqz fc log(2;) + 0{z fc), 


( 18 ) 


with n = kp, we obtain : 




k 


WP-1 \ r 

^^^Lvol(0)j J \t\P-^g{t,Xo)dt. 

' ' R 


c) If n > k and n/k ^ N (integrable singularity) we obtain the same result as 
in a) but with the modified distributions : 


tf ^,g{t,xo)) (-^wl \Lvo\\ + (t!l ^,g{t,xo)) (-^wl ^ Lvol), 


dw^- 


dw^ 


where the derivatives w.r.t. w are normalized distributional derivatives. 

The meaning of normalized derivative it that one choose the normalization 


, dw' 


-w± ‘‘,f{w)\:=(w±‘‘,f{w)j, 


(19) 


for all / G with / = 0 in a neighborhood of the origin. The distributional 
bracket involving Lvol is detailed in section 3. Results c) and b) for p > 2 
are not intuitive and are certainly difficult to be reached without geometry. 
In particular, for applications to oscillatory integrals (see below) one has to 
work in the dual since both Fourier transforms w.r.t. t in c) and b) are distri¬ 
butional. In c), the n-th derivative is arbitrary and the result is the same for 
any normalized derivative of order greater than E(n/A;). 

Remark 6 Results a) and b) for p = 1 are certainly interesting for spectral 
analysis since these contributions are bigger than Ir{z) = 0{z~^). As in The¬ 
orem 3, non-integrable singularities have a dominant contribution. Hence, the 
leading term of I{z) is always an invariant. 

Application to oscillatory integrals. 

A typical application of Theorems 3&5 can be the asymptotic expansion of 
distributional traces of quantum propagators. Hence, it is interesting to remark 
that our results can be extended to asymptotic integrals : 

I{z) = J G{z, zf{x), x)dx, z ^+ 00 , 

X 
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if G admits an asymptotic expansion with a priori estimates, i.e. : 

i 

G{z,t,x) = gj{t,x) + Ri{z,t,x), 

j=0 

ykeN*: \\Rk{z,t,x)\\mR>,x) = £ > 0, 

where {aj)j is a strictly increasing sequence and R controlled by uniform 
estimates. Similarly, we can consider expansions in term of z~°‘^ log{z)"^. This 
notion of graduation w.r.t. allows to apply our results but, to simplify the 
exposition, in this work we will just consider the case of an integral of a symbol 
g(t,x). We can treat degenerate oscillatory integrals : 

0(z) = J z)dtdx, z ^+oo, (20) 

RxX 


providing that / satishes the conditions of Theorem 3 or 5. An important 
application in quantum mechanics is the case X = T*W^ where, after some 
technical modihcations, the localized (distributional) trace of h-pseudors : 

TrUhiAh - E) := Tr J e C'o“(M), E eR, 

R 


can be written as a locally hnite sum of oscillatory integrals : 


J e^^^^^'y''^^Ayx)-tE)b^hR,y,ri)dtdydri, (21) 

RxT*R" 

where h{h, •) ~ I] h~%k satishes a priori estimates as above and S is the local 
generating function of the group of diffeomorphism of the principal symbol of 
Ah. Here z = h~^ is the parameter and, after a discussion based on classical 
mechanics, Eq. (21) can be reformulated as in Eq. (20) where © is the energy 
surface of level E. For more details, we refer to [2,4,5]. 

We recall now basics on homogeneous transformations, some of them will also 
be used below. The Melin transform of a function h is dehned as : 

OO 

M[h](0 = J h{t)t^-^dt. (22) 

0 

Generally M[h]{^) is analytic, for example, in the strip 3?(,^) E]a,b[ with : 

OO OO 

a = inf{ / \h{t)\R~^dt < cx)}, b = sup{ [ \h{t)\R~^dt < cx)}, 
xeR J a;e]R J 
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and when M[h]{c + iy) G L^(M, dy) we have the inversion formula : 

C+200 

= ^ (23) 

c —200 

By elementary changes of path we obtain that : 

= e^“^r(0, Re(0 e]o,i[. 

We write 0(z) = 0+(z) + 0-(z), where : 

0±(z) = J e^^^^^^^a{t,x)dtdx. 

{±tf{x)>0} 

Melin’s inversion formula leads to the distributional formulation : 

C+200 

0±(2:) = — J J \tf{x)\~^a{t,x)dtdxd^, 

c-ioo {±tf{x)>0} 

where 0 < c < <5 for h > 0 depending on /. If no elementary method works, 
one can construct explicit meromorphic extensions of distributions : 

J {tf{x))~^a{t,x)dxdt, 

{tf{x)>0} 

(R/“(0,«) = J \tf{x)\~^a{t,x)dxdt, 

{tf{x)<0} 

and use Cauchy’s residue formula to get asymptotic w.r.t. z ^ oo. This 
method is described in [8] when / is monomial, which is a generic situa¬ 
tion when / is analytic on supp(a) by Hironaka’s theorem of resolution of 
singularities [6]. Also in many cases it is possible to reach this setting by a 
non-analytic change of coordinates. This approach is very interesting but a 
great disadvantage is that there is 4 domains, a sequence of poles (of order 
3) which does not contribute at all in the expansion and the complex factors 
exp(±i7r.^/2)r(,^) lead to long calculations in presence of multiple poles. 

In fact we will reach Theorem 3&5 more directly and by a method which 
avoids Fourier analysis until the last step of the proof. 


2 Case of © regular and non-degenerate critical points. 

If / is regular on ©, then © is a compact and smooth submanifold. 
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Proposition 1 Under the previous assumptions and if & is a regular surfaee, 
I{z) admits a full asymptotie expansion in powers of z~^ with : 

= “ / / 9 (t,x)dtde{x) + 0{z~‘^), (24) 

x£6 K 

where d^ is the invariant surfaee measure of &. The same result holds for the 
integral Ir{z) with insertion of the eut-off in the integral. 

By the implicit functions theorem and compactness, we pick some coordinates 
y, dehned near ©, and a hnite partition of unity Qj covering ©, such that / 
is diffeomorphic to the coordinate yi in supp(f2j). By Lemma 2 we obtain ; 

/w = y/,w + o(2-“), 

re) = J j 

R 2 R ’^-1 


where y* and \ Jy\ are respectively the pullback by y and the multiplication 
by the standard Jacobian. Hence, we have to estimate an oscillatory integral 
with quadratic phase. The stationary phase Lemma, see e.g. Lemma 7.7.3 of 
[7] vol.l, provides a full-asymptotic expansion : 


/ e'^^'^^y*{VLj{x)g{t,x)\Jy\)dtdyi = 

R2 


N-1 

1=0 


(25) 


Since sign(t|/i) = 0, the expansion is real with leading term : 

Stt /* 

^i(2:) = — j iy*i^jix)git,x)\Jy\)){0,0,y2,...,yn)dy2--dyn + O{z~‘^). (26) 

R"-l 

Geometrically, since & is locally given by yi = 0, this corresponds to integra¬ 
tion on © w.r.t. the canonical n — 1 dimensional measure de. By summation 
over the indices j we obtain the integral on the whole surface, i.e. : 

I{z) ~ — y 9 { 0 ,x)de{x) = - J J 9 {t,x)dtde{x). (27) 

x£& x£6 R 


Remark 7 Eq. (27) is a particular (flat) case of stationary phase formulas 
with a smooth compact manifold of critical points and non-degenerate trans¬ 
verse Hessian, d© is the Liouville-measure of classical mechanics or Guelfand- 
Leray-measure in theory of singularities. The oscillatory representation of 
delta-Dirac distributions, by mean of Schwartz kernels, provides a very natural 
dehnition of this object. 



Now, assume that / admits a single non-degenerate critical point on &. As 
usually, non-degenerate means that Q{^) = | (^, (Pf{xo)^) is a non-degenerate 
quadratic form. After perhaps a reduction of the cut off 0, the Morse Lemma, 
with coordinates w, provides a simpler problem : 

F{z) = j F^^^^'^'’w*{g{t,x)Q{x)\Jw\{x))dw. (28) 

KxR" 


This kind of asymptotic problem was precisely studied in [2], The result is ; 
Proposition 2 If Q has at least one even index of inertia then : 

OO 

(29) 

j-o 

But, if Q has both indiees of inertia odd then we have : 

OO 

F{z) ~ ^ ^ Cj^mZ~^-nog{z)'^. (30) 

jr=0 m=0,l 

The method used in section 3.2 allows to find this result. Our proof is not 
simpler and that’s why we refer to Proposition 3.4 and Theorem 3.5 of [2] for 
a detailed proof and discussion of the coefficients. This discussion is necessary 
because the exponential distributions of Eq. (28) possess spherical symmetries. 


3 Proof of the main results. 


3.1 Critieal points attaehed to a local extremum. 


We start by the case where / admits a local minimum xq on the fiber ©. We 
split-up our integral I{z) via Eqs. (9,10). A can be treated as in section 2 and 
from now we concentrate our attention on A. With the extremum condition, Xq 
is isolated on © and ©nsupp(0) = {xq} for supp(0) small enough.To simplify 
notations we identify xq with the origin and we recall that we consider : 

fix) = fk{x)+ 0{\\x\\’"^^), (31) 

where fk is homogeneous of degree k and dehnite positive. Using polar coor¬ 
dinates, the Taylor formula shows that : 

/(r9) = r‘(/t(9) + fi(r,9)), (32) 
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with i?(0, 6^) = 0 and where we note again fk{0) the restriction of fk to S” 

If supp(0) C -8(0, To) is chosen small enough we have : 

fk{9) + R{r,9) ^ 0, e Vr e [0,ro[. 

Hence, with the homogenous coordinates v = {u, 9) where : 

u{r-,9) = r{fk{9) + R{r-,9))^, (33) 

we can express our integral as : 

Is{z) = ( G{zu^,u)du. (34) 

R+ 

The new symbol G is obtained by pullback and integration : 

G{t,u)= j v*{g{t,r9)Q{r9)r^-\j{v)\)d9. (35) 

§n-l 

The existence of the asymptotic expansion is a consequence of : 

Lemma 8 For a in C“(M x M+), the following asymptotic expansion holds : 

J{z) = / a{zu^,u)du 'dj{a), z ^ oo, (36) 

i>o 


where the dj are universal distributions given by : 


(j. — ^ ^ 


® E («))• 


(37) 


Proof. The crucial point is that is increasing on M+. We have : 

OO 

J{z) = z~k f a{u'^,—)du. 

J Zk 
0 

A Taylor expansion w.r.t. the second argument at the origins gives : 

^ 1 i ,d^a 


a(M^,E)=EE + z Rn+i{u,z), (38) 

zk i^qU ow 

where Rn+i{u, z) is integrable w.r.t. u with norm uniformly bounded in 
By a new change of variable we obtain : 


Y/ \ 1 V—^ 1 _ 1+^ r d a . . i+i—k N+i 

= ' dr + C>(^ ). 

1=0 ' n 


(39) 
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These coefficients are well defined since |r e LLm for all I eN. m 

Remark 9 This expansion holds also for a pullback by —u^ if we replace r" 
by r". This allows to treat the case of a local maximum in Theorem 3. Also 
for an application to oscillatory integrals we obtain a nice formulation via the 
Fourier transform of the distributions r^, which avoids any ’’regularization”. 


Proof of Theorem 3. We treat the case of a local minimum. We apply Lemma 
8 to the integral of Eq. (34) to prove the existence of the asymptotic expan¬ 
sion and it remains to express invariantly the leading term. With the polar 
coordinates G vanishes up to the order n — 1. Consequently, we have : 


Liz) = 


[n 


1 )! 


'i — k 
k 




Starting form Eq. (35) and since : 

\jv\{o,e) = \Me)\-i, 


by elementary manipulations on delta-Dirac distributions, we obtain that : 

L{z) = J \fk{0)\~'^d9 + O{z~'^). (40) 

gn-l 

n—k n—k 

Finally, for a local maximum we replace the distributions by tL . ■ 

Remark 10 The leading term is invariantly defined since Eq. (40) involves 
the evaluation g{t, Zq) and fk is an invariant since the [k — l)-jet in zq is fiat. 


On the integral on the sphere. 

The integrals on of Eq. (40) can be reformulated. For example in dimen¬ 
sion 2 these are elliptic integrals. In higher dimensions we define : 

Hfk) = I (41) 

since fk is homogeneous we obtain : 

OO 

I{fk)= J = J e-^'u^-^du J \fk{e)\-^d9. 

R+x§"-i 0 §"-i 

Hence our integral is given by : 

§n-l ^ I > 
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Iifk) can be computed, by elementary methods, as a product of gamma factors 
or hypergeometric functions. We give an elementary example in section 4. 

3.2 Case of non-extremum critical points. 

This problem is more complicated and we recall that we just treat here sin¬ 
gularities given by {Hi). As before, we just have to study Is{z) and, in local 
coordinates, we identify the critical point with the origin. We dehne C{fk) as 
the trace on the unit sphere of the conical set of the zeros of fk, i.e. : 

C{fk) = {9e / Me) = 0} = n {x e RV fk{x) = o}. (43) 

With {Hi), C{fk) is a compact smooth submanifold of of codimension 
1. To perform a blow-up of the singularity we use polar coordinates (r, 6) and 
the next lemma gives a resolution of the singularity w.r.t. C{fk). 

Lemma 11 In a micro-local neighborhood of the origin there exists local co¬ 
ordinates y, on the blow-up of the critical point, such that : 

/ 

y\, in all directions where fk{e) > 0, 
f{x) — \ —i/i, in all directions where fk{e) < 0, 
yM, locally near C{fk). 

Proof. By Taylor, there exists R continuous in r = 0 such that : 

f{x)c^f{re)=r\Me) + R{r,e)). (44) 

If 6*0 ^ C{fk) and 9 is close to 6 *o we simply choose : 

{y2,...,yn){r,9) = (0i,...,0„_i), 
yi{r,9) = r\fk{9) + R{r,9)\P 

In these coordinates the phase becomes if fk{9o) is positive (resp. —yi for 
a negative value) and the Jacobian satishes \.Jy\{t),9) = |/fc( 6 *)|fc ^ 0 locally. 
Now, let 6*0 G C{fk). Up to a permutation, we can suppose that i9ei/fc(6'o) 7 ^ 0. 
We accordingly choose the new local coordinates : 

{yuy3,-,yn){r,9) = (r, 02 ,--,^n-l), 
y2{r,9) = fk{9) + R{r, 9). 

Since we have |J|/|(0,6*o) = \d0j^fk{9o)\ ^ 0, lemma follows. ■ 


12 



Remark 12 These coordinates are admissible and this leads to an adapted 
system of charts via a partition of unity on introduced below. Also the 
condition {Hi) insures the existence of a canonical measure on C{fk)- 

To use Lemma 11 we introduce an adapted partition of unity on We 

pick cut-off functions e 0 < < 1, X) = 1 in a tubular 

neighborhood of C{fk), with supports chosen so that normal forms of Lemma 
11 exist, for r small enough, in a conic neighborhood of supp('Lj). By com¬ 
pactness this set of functions can be chosen hnite and we obtain a partition 
of unity on by adding Tq = 1 — X fo our family. The support of Tq is 
not connected and we dehne Tq , with fk{0) > 0 on supp('l/g), and similarly 
we dehne Tg where fk < 0, so that Tg = Tg -|- Tg . If we accordingly split up 
Is{z) we obtain : 

= J '^Q{6)g{zf{r6),rd)Q{r6)H~^drd6 

RxR+XS2n-l 

= / Go{^zyi,yi)dyi, 

R+ 


respectively for the directions where fk{0) > 0 and fk{0) < 0, also : 

I^/{z)= J '^j{e)g{zf{re),re)Q{re)H-^drde 

RxR+x§"-i 

= J Gj{zy\y2,yi,y2)dyidy2, 

R+ xR 


for the set G{fk)- The new symbols are respectively given by : 

aiit.m) = I v‘{'l'i{l>)g{t,re)B(r0)r’‘-'\.Ji,\)dy,...dy„, (45) 

Gj(t,yuy2) = J V'‘{'i^A0)9(*<r0)O(re)r’‘~'\Jy\)dy3...dyn- (46) 

Hence, the singular part of our integral can be written as a hnite sum : 

I.(z) = I-(z) + I*(z) + x: If’iz), (47) 

j 

where each term of the r.h.s. will be treated by elementary methods. Note 
that I~ {z) and {z) can be treated as in the previous section. 

Remark 13 Since yi{r, 6) = r, our new symbols satisfy Gj{t, yi, 1 / 2 ) = G{yi~^), 
near yi = 0. Since the asymptotic expansion involves delta-Dirac distributions 
w.r.t yi, cf. Lemma 14 below, the dimension will cause a shift in the expansion. 
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For a G (^“(M x M+ x M), we define the family of elementary fiber integrals : 

OO 

InM = /(/ a{zy\y2,yuy2)dy2)yi~^dyi. (48) 

0 R 

Lemma 14 There exists a sequenee of distributions {Dj^p) such that : 

dn,k{^) ~ H Dj^p{a)z~i log( 2 :)P, as z ^ oo, (49) 

p= 0 ,l jEN, j^n 

where the logarithms only occur when {j/k) is an integer. As concerns the 
leading term, if{n/k) ^ N* we obtain : 

In,k{z) = z~^d{a) + 0{z~^log{z)), (50) 

with : 

OO OO 

d{a) = Cn,k J J 0 , 1 / 2 ) + dp^a{-t, 0, - 1 / 2 )) dy 2 dt. 

0 0 

But when n/k = p eN* , we have : 

dn,k{z) = ^z~'^\og{z) J \tf~^d^f^a{t,0,0)dt + O{z~'^). 

R 

7X I 1 

Remark 15 The remainder of Eq. (50) can be optimized to 0{z~~iz ~) when 
(n + l)//c is not an integer, as shows the proof below. 

Proof. By a standard density argument we can assume that the amplitude is 
of the form a{s, yi, 1 / 2 ) = f{s)b{yi, y 2 ). The justification is that our coefficients 
below are computed by continuous linear functionals, i.e. distributions. We 
define the Melin transforms of / as : 

OO 

M ±(0 = J s^-^f{±s)ds. (51) 

0 

We split-up In,k as /+ and /_ by separating integrations 1/2 > 0 and 1/2 < 0. 
Via Melin’s inversion formula, we accordingly obtain : 

I+{z) = ^ J M+{f)z-^ J {yiy2)~^b{yi,y2)yi~^dyidy2df, (52) 

7 R^ 

where 7 = c + iM and 0 < c < k~^. Similarly we have : 

I-{z) = ^ J M_{^)z~^ J {yiy^y^biyi, -y2)yi~^dyidy2df, (53) 

7 R2 
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The existence of a full asymptotic expansion is a direct consequence of : 

Lemma 16 The family of distributions i—>• on initially 

defined in the domain < k~^ is meromorphie on C with poles : = j/k, 

j G N*. These poles are of order 2 when ^ Id* and of order 1 otherwise. 


Proof. We form the Bernstein-Sato polynomial bk attached to our problem : 

T{y2y’iy~^ := Qy^Qyk (^2l/i)^~^ = bk{0{ty2yi)~^, 

k 

j=i 

If < /c“\ {k + l)-integrations by parts yield : 


/ (1/1I/2) ^f{yi,y2)dyidy2 

R 2 




J ( 1 / 11 / 2 )^ HTf)iyi,y 2 )dyidy 2 . 


Now the integral in the r.h.s. is analytic in < 1 + k After m iterations 
the poles, with their orders, can be read off the rational functions : 

m -I 

This gives the result since m can be chosen arbitrary large. ■ 

Remark 17 Starting from Eqs. (52,53) it is possible to obtain an asymptotic 
expansion as O’*". By shifting the path 7 to the left, only the poles of the 
Melin transforms contribute and we obtain an expansion in power of 2 :. This 
can be predicted by a Taylor expansion of the associated oscillatory integral 
but the improvement is that one can establish sharp integral remainders. 

Hence, the following functions are meromorphie on C : 

0^(0 = J {yiy2)~^b{yu±y2)dyidy2. (55) 

R 2 


A classical result, see e.g. [1], is that M±(c + ix) G 5(Ma,) when c ^ —N. If 
we shift the path of integration 7 to the right in our integral representation, 
Cauchy’s residue method provides the asymptotic expansion. In fact for any 
d > c, outside of the poles, we have : 

I+{z)= J 

c+iR 

= ^ res(W^M+g)(^j-fc) + J W«M+(O0+(Oc?^- 
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Since d is not a pole the last integral can be estimated via : 


1 J z-^M^{05^{0d^\<C{f,b)z-‘^ = O{z-^[ 

d+iR 

), (56) 

where, for each d, the constant C involves the L^-norm of a finite number 
derivatives of b. This will indeed lead to an asymptotic expansion with precise 
remainders. Applying this method to I+{z) and I-{z) we obtain the existence 
of a full asymptotic expansion of the form : 

Hz)~ E E c„2-iiog(zr. 

p=0,l jGN* 

(57) 

Moreover, by Lemma 16, these logarithms only occur when j/k is integer. 

Now, we compute the leading term of this expansion for our 
lem. To avoid unnecessary discussions and calculations below, 
we can commute the polynomial weight via : 

particular prob- 
we remark that 


(58) 

6(0 = + "-!)■ 

(59) 


t=i 

By iteration, we obtain that the poles are the rational numbers : 


1 Tl — 1 

^p,j,k,n = P+ -^, j e[l,...,k], pe N. 

For all a — /3 > 1, /3 G N*, we have : 

OO 

J d^{ry{r)dr) = 0, V/ G 

0 

If we apply this to the integral w.r.t. yi, there is no contribution before : 

eo = I ( 60 ) 

To compute the first effective residue we must distinguish out the case where .^o 
is an integer or not. The optimal number of iterations to reach is E(n/ k) + l 
but, by analytic continuation, any number bigger than this one is acceptable. 
A fortiori we can use n iterations and our starting point will be ; 

f (yiy2r-^v'!-'T"b{m,ih)dmdy,, (61) 


with : 


n—1 

»«(o=n 


1=0 


1 


(62) 
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Case of ,^0 simple pole. 

In this case our residue is simply given by : 

J {y^y2)'"~'^y'l~^T^b{yi,y2)dyidy2, (63) 

R 2 

with : 

lim(-ir''=+‘'K-?)»n(0. (64) 

k ^ 

In particular we can compute the integral w.r.t. yi via : 

OO 

/= (-!)*"((;« - l)!S“M 0 ,t/ 2 ). 

0 

A similar result holds for /_ and we obtain : 

OO 

n n P n—H 

I+{z) = z-kCn,kM+{-) J y 2 " (<9” &)( 0 , y 2 )dy 2 + Ri{z), (65) 

0 

OO 

71 r ^ 

I-{z) = z~'^Cn,kM-{-) J y 2 {dy^b){0, -y 2 )dy 2 + R 2 {z). ( 66 ) 

0 

Here Cn,k is the canonical constant : 

Cn,k = T n n- (67) 

i=l d k 

7X I 1 

Also, according to the analysis above, each remainder is of order 0{z ~) if 

77-^h 1 

(n + 1 )//;; ^ N and 0{z ~ log(; 2 )) otherwise. 

Case of ,^0 double pole. 

If h is meromorphic with a pole of order 2 in ,^o we have : 

res(h)(^o) = 

Applying this principle to our residue we obtain, via Leibnitz’s rule, that : 

I+{z) = B \og{z)z~'k + 0{z~'^). (68) 

We can compute the distribution B as before and we hnd : 

1 °° 

I 71 C rt — R 

B = --D„,,M+(-) J yrHd^M0,y2)dy2, 

0 
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Since p = n//c is an integer, by integration by parts we obtain : 

OO 

/!/r’'(aj6)(0,j/.)<ij/2 = (-l)’-''+‘(n-p)!aj;‘M0,0). (69) 

0 

Since a similar result holds for we obtain the desired result by gathering 
all the constants and summation. Finally, we can extend our formulas since 
all coefficients in the expansion are of the form : 

{Ti,f»b) = {t{, /) {Ti,b) , Tl, € V. 

By linearity and continuity, the result holds for a symbol a(t,yi,y 2 ). ■ 

Taking Remark 13 into account, to avoid unnecessary calculations we define : 

Go{t, yi) = yi)^ (70) 

Gj{t,yi,y2) =yi~^Gj{t,yi,y2). (71) 

Directions where fk{0) ^ 0. 

By Lemma 8, the hrst non-zero coefficient, obtained for / = n — 1, is : 

n. _ n 

n — k / 1 \ \ 2' k r f^—k ^ 

t~ 0 5^”“ \G^{t,yi)j = J t~G^{t,0)dt. 

R 

By construction, we have : 

G„+(i,o)= I g(t,(i)<ii(e)\M 0 )\-ide. 

§n-l 

A similar computation gives the contribution of supp(\kQ ), and we obtain : 
i:{z)^z-i(t2,9{t,0))l I n(Wzm-iM+0(z-"-i^), (72) 

§n-l 

/.-(*-) = *--* (iW,g(i,0)) 1 / %mMe)\-’^de+0(z-"-^). (73) 

§n-l 

Microlocal contribution of G{fk). 

Here we examine the contribution of terms Ig’^{z). According to the analysis 
above, we must distinguish out the case k divides n. 

a) Case of k > n, non-integrable singularity on &. 

Here n/k g]0,1[, so that the singularity on the blow-up is integrable. Via 
Lemma 14, the contribution of Ig’^{z) is given by ; 

^z~'^ J |t|^“^|p2r^ (Gj{t, 0 ,y 2 ) + Gj{-t,0, -P2)) dtdy 2 + 0{z~^\og{z)). 

Rl 


Z k 1 

k (n — 1)! 
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Reminding that y 2 {t, 0, 6) = /fc(6'), we obtain : 


y I 1 / 2 I fcGj(t,0,|/2)c^l/2 = fi'(t,0) J \fki0)\ ^^ji9)de. 

*+ {/fc( 0 )>O} 


Since \1/^, \I/j is a partition of unity on S” by summation of all local con¬ 
tributions Is{z) is asymptotically equivalent to : 


_ 71 

2 ; k 


k 


{Ux \^(t,0) 


V 


J \fk{9)\ kde+(t!i \g{t,0) 

{/fc>0} 


/ IfkieT'^de 

{/fc< 0 } 


Note that none of these coefficients are equal unless g or are symmetric. ■ 
b) Case of p = n/k integer. 

Here the contribution of is dominant and we obtain : 

Ipiz) ~tlog(z)z-f I \t\’’-'dlpG,{t,0,(l)dt+ 0{z-n. 

R 


Unless p = 1, there is no way to take the limit directly, and the geometric 
properties are still hidden in the Jacobian. But we will reach the result by the 
Schwartz kernel technic. Clearly, it is enough to evaluate our derivative and 
to integrate w.r.t. t. With s = (si, S 2 ) G M^, we write the evaluation as : 

J ed^'^y^'^^^\is2y~^Gj{t,yi,y2)dyidy2ds. 

Here we have used an oscillatory Schwartz kernel for Sy.^ ® This integral 
representation allows to inverse our diffeomorphism to obtain : 


d^yf^Gy(t,0,0) 


1 


ed^dr,y2(r,m^is^y-itP-ig^t^r9)^.^g^dtdrd9ds. 


Extending the integrand by 0 for r < 0, the normalized integral w.r.t. (r, Si) 
provides dr- By construction 7 / 2 ( 0 , 6*) = fk(9), hence : 


^y2 = 9(i,0) 





(74) 


This Fourier integral makes sense with ^ compact. We recall the density : 


Jj(w) = I ^,{e)dLU0), (75) 

{fkW=w} 

where is the density induced by the Leray-form dLjj, : dfk/\dLf^{6) = dO. 
Note that all these objects can be constructed by mean of local coordinates 
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under the only condition that snpp(\l/j) is small enongh near C{fk). Since fk 
is continnous on each Jj{w) is compactly supported and smooth near 

the origin. The sum over all the gives the geometric contribution : 



EM w)dwdu 
j 


dP-^Lvol 

dwP~^ 


( 0 ). 


(76) 


By integration w.r.t. t we obtain the result for a double pole in general position. 


Remark 18 The case n = k is directly accessible. In this case we obtain : 
I^(z) = ^i^^LVol(O) Jg{t,0)dt + O{z-^). 

R 


Here LVol(O) is the Liouville volnme of C{fk)- Note that Ig dominates J^. 
c) k < n and simple pole, integrable singularity on &. 

Finally, we treat the case of a simple pole with a non-integrable singnlarity on 
For the positive part, the distributional coefficients are : 

OO OO 

{t^+,9) = Cn,k J J t'^~^yT^{dy^Gj){t,0,y2)dy2dt. (77) 

0 0 

Clearly, we can use the same globalization technic as above. The snm of all 
completed with the main term of J+, provides : 

OO OO f)nr ( \ 

{T+,9) = Cn,k J J (78) 

0 0 

where the integral w.r.t. w is absolutely convergent since the measure L{w) is 
compactly supported. A similar result holds for the directions where fk{0) < 0 
and we obtain the resnlt stated in Theorem 5. ■ 

Now we detail the construction of the distributional bracket of the part c) of 
Theorem 5. Let be y G C^, 0 < y < 1 on R, chosen such that x = 1 near the 
origin and x(m) = 0 for |u| > e. If e > 0 is small enough, Lvol is smooth on 
] — e,e[. We write the geometric contribntion as : 

(T, Lvol) = (T, xLvol) + (T, (1 - x)Lvol) 


Away from the origin, e.g. for w > 0, we obtain directly : 


- x)(m)Lvo1(u) 


\ I (i-xKhimhwMM. 

{fkid)>0} 

(79) 
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Note that Eq. (67) for Cn,k justifies the normalization of Eq. (19). On supp(x), 
we use the local regularity of Lvol(u) near u = 0 and integrations by parts to 
conclude. Finally, to obtain a totally rigorous treatment, we remark that since 
C{fk) is compact and fk is continuous we can choose our partition of unity 
such that Z) = 1 for \fkiO)\ < 2e. ■ 

Remark 19 The key point here is that we can put in duality the distributions 
and Lvol(u) since their singular supports are disjoints. 


4 Elementary examples. 


A family of extremum. If p is an even integer and a = (oi, ...an) G (M!j_)"', 
we consider a n dimensional integral obtained by integration in the hber of : 

n 

fa,p{^) = ( 80 ) 

i=i 

For example ; 

n "1 ^ 

e-^f^,p(^)dx = n(/ = (2r(l + -))^z-p H S ^ 

K" i=i R P i=i 

with g{t,x) = symbol of order 0 on R”. Since we have here a local mini¬ 
mum, the distributional factor is simply : 


/ n — p , 

t p e~^dt = r{n/p). (81) 

0 


The Jacobian of the standard polar coordinates gives the value of the gener¬ 
alized elliptic integrals : 


gn-1 


1 (2r(i + P)) 

p T{n/p) 



(82) 


Note that Eq. (82) can be analytically continued with a phase Z « > 0. 

Conical singularities. In 2 dimensions let be : 


g{t,u) 


e-IMP 

(l+^‘ 
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With f{x,y) = we have © = {(x, x)} U {(—x, x)}. © is not compact 

but the exponential decrease of g will compensate. Accordingly let be : 



1 + z‘^{x‘^ — y'^Y 


dxdy. 


(83) 


Since n/k = 1, we expect a logarithm for the leading term. Since : 

27r 


dO 


271 


0 


1 + a cos^(26') + a 


, Va > —1, 


by passage in polar coordinates we obtain : 


OO 

I{z) = 271J e~ 


rdr 


du 


vT+ 


2r4: 


= 7T e 


z^r 


vl + 


z'^u'^ 


TT- 


.log(^) 


Always in dimension 2, a singularity of degree 4 will generally not give a 
logarithmic leading term. For example, we have : 


i{z) = 


3-(a;2+j/2) 


1 + z‘^{x^ — y^^Y 


dxdy = 71 e~ 


du 


7T 


du 


Vi + / vTT' 


u^ 
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